Hunt et al. Neural Systems & Circuits 2011, 1:3 
http://www.neuralsystemsandcircuits.eom/content/1/1/3 



7g\ Neural Systems & Circuits 



RESEARCH Open Access 



Statistical structure of lateral connections in the 
primary visual cortex 

Jonathan J Hunt 1 , William H Bosking 2 , Geoffrey J Goodhill 1 * 



Abstract 

Background: The statistical structure of the visual world offers many useful clues for understanding how biological 
visual systems may understand natural scenes. One particularly important early process in visual object recognition 
is that of grouping together edges which belong to the same contour. The layout of edges in natural scenes have 
strong statistical structure. One such statistical property is that edges tend to lie on a common circle, and this 
'co-circularity' can predict human performance at contour grouping. We therefore tested the hypothesis that long- 
range excitatory lateral connections in the primary visual cortex, which are believed to be involved in contour 
grouping, display a similar co-circular structure. 

Results: By analyzing data from tree shrews, where information on both lateral connectivity and the overall 
structure of the orientation map was available, we found a surprising diversity in the relevant statistical structure of 
the connections. In particular, the extent to which co-circularity was displayed varied significantly. 

Conclusions: Overall, these data suggest the intriguing possibility that VI may contain both co-circular and anti- 
cocircular connections. 



Background 

Grouping edges which belong to the same object or 
contour is a vital part of object recognition. Biological 
vision systems excel at this task, yet it is still extremely 
challenging for artificial vision systems [1]. In humans, 
contour detection develops well after birth [2], suggest- 
ing that learning from visual experience plays an impor- 
tant role. 

The layout of edges in natural scenes has strong statis- 
tical structure. One well-known property is co-linearity. 
Co-linearity is important in contour grouping and tex- 
ture discrimination tasks [3-5]. A more recently discov- 
ered generalization of co-linearity displayed by natural 
scenes is co-circularity: edges tend to be tangent to a 
common circle more often than would be expected by 
chance [6]. Furthermore, the degree of co-circularity in 
a contour can be used to predict human contour detec- 
tion performance [7]. This raises the question of what 
biological substrate underlies this effect of co-circularity 
on performance. 
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An obvious candidate is the pattern of excitatory lat- 
eral connections in the primary visual cortex (VI). 
These connections are primarily found in layer 2/3 and 
are longer range than the inhibitory connections found 
predominately in other layers [8]. In humans these con- 
nections develop after birth, and their development 
coincides with improvement in contour detection 
[2,9,10]. Furthermore they have been shown to play a 
role in co-linear facilitation [11-14]. 

It has been previously demonstrated that long-range 
lateral connections connect preferentially to patches of 
similar orientation and primarily along the axis of orien- 
tation in VI [15-18] (although others have claimed this 
preference is weak [8]). This is consistent with their 
role in co-linear facilitation. Although previous work 
has demonstrated that the variance in some of these 
data is not adequately explained by co-linearity alone 
[19], whether there is a more general tendency towards 
co-circularity in these connections has not been 
investigated. 

Here we reanalyzed previously published functional ana- 
tomical data regarding lateral connections in tree-shrew 
[17], using a noise-resistant measure of co-circularity 
that we recently introduced to study the structure of 
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orientation maps [20]. Surprisingly, we found a large varia- 
tion in the statistics of lateral connections between ani- 
mals. One unifying explanation for both this and previous 
results is that lateral connections in VI are connected in a 
variety of ways, both co-circularly and anti-cocircularly. 
We suggest reasons why this may be a desirable arrange- 
ment for analyzing the structure of visual scenes. 

Methods 

We reanalyzed data from four tree shrews previously 
published by Bosking et al. [17]. The experimental 
methods used are described in detail in that paper and 
we provide only a brief overview here. Tree shrew orien- 
tation preference maps were obtained using optical ima- 
ging. Additionally, 540 nm light was used to map 
surface blood vessels used for alignment. Biocytin was 
then injected into a specific site in VI and the animal 
was sacrificed 16 hours later. Slices of VI were imaged 
to locate the biocytin bouton and the surface blood ves- 
sels. The blood vessel information was then used to 
align the orientation preference maps with the bouton 
images giving overlaid information on the underlying 
connectivity from the injection site on the animal. The 
original experiment used a total of ten cases, however, 
we were only able to recover the data for four cases. 

Topography 

In order to accurately quantify the co-circularity present 
in these results it is important to know the underlying 
topography. Fortunately, the topography of the tree- 
shrew is consistent between individuals and well charac- 
terised [21]. To align the maps (and overlaid boutons) 
with visual space they were first rotated so the V1/V2 
border was vertical, then flipped along the vertical axis 
to give a right-handed co-ordinate system (shown sche- 
matically in Figure 1). The tree-shrew VI has a com- 
pressed representation of the ipsilateral visual field near 
the V1/V2 border with a very different magnification 
factor from the rest of VI. We therefore drew a line 
representing the edge of the contralateral visual field 
(the vertical meridian) and eliminated the ipsilateral por- 
tion of the map. The map was then rotated slightly to 
ensure the vertical meridian was represented vertically. 

This gave us an orientation preference map of VI with 
an aligned injection site and bouton sites which indicate 
the lateral connections originating from the injection site. 
The rotating and flipping ensured that points on the map 
corresponded to points in visual space within a scaling fac- 
tor. After truncation of the ipsilateral hemisphere the mag- 
nification factor was approximately constant for the region 
of VI imaged, and cortical distance and direction could be 
used as a reliable proxy for visual field distance and direc- 
tion. The aligned maps with overlaid bouton and injection 
sites are shown in Figure 2. 



For the original analysis [17] electrophysiological 
recordings taken during the injection were used to 
determine orientation. We determined the orientation 
of the injection site using the orientation value of the 
underlying pixels that were measured using optical 
imaging. This differed from the electrophysiological 
recordings by an average of 15° (range 0-23°) from 
the optical orientation value (averaged across the 
animals). 

Characterizing the lateral connection statistics 

To test for the presence of co-circularity in the lateral 
connections we adapted a noise-resistant measure of 
co-circularity introduced in Hunt et al. [20]. For each 
pair of injection/bouton sites on the orientation prefer- 
ence maps we defined the following (also shown dia- 
grammatically in Figure 3A): 

9 i orientation preference of injection site i; 

9 b orientation preference of bouton site b; 

cortical distance between injection site 

r 

and bouton site; 

angle of the line from injection site i 
to bouton site b; 

relative orientation preference of the bouton b 

a =6 h -8, (mod 180°) 

relative to injection site i; 

angular position of the bouton site b 
/J=^-e,-(modl80°) f ■ u • 

relative to the orientation prelerence 

of the injection site i. 

The orientation preferences 6 t and 6 b are defined in 
real space where the stimulus is presented, relative to 
zero degrees along a horizontal #-axis. They are always 
taken on the j-positive side of the horizontal and thus 
range from 0° to 180°. The relative position of the corti- 
cal points 4> is defined relative to the horizontal axis of 
the retinotopic map and is calculated modulo 180°. 
Except where noted, calculations of orientations and 
angular positions were done using modulo 180° 
arithmetic. 

The wavelength A of each map was calculated as the 
mean of the Fourier spectrum averaged over all direc- 
tions as previously described [22]. As there is variation 
in column spacing between individuals [23], and the 
imaging of each map may be at different magnifications, 
the use of A as the unit of length aims to ensure that 
the quantification of reduced symmetry is not affected 
by overall changes in map scaling. The wavelength gives 
a rough measure of the size of iso-orientation patches in 
the map. 

We calculated a probability distribution for the con- 
nection statistics for each animal p(a , [}, r). However, if 
(f> is not uniformly distributed this can introduce spur- 
ious correlations. We corrected for this by calculating a 
four-dimensional probability p(a , ji, r, cf>) and then 
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Figure 1 Topography matching of tree shrew maps. (A) shows a schematic of the topography of the tree shrew. The area shaded grey is a 
compressed topography of the ipsilateral hemisphere and is excluded from analysis due to its very different magnification factor. To facilitate 
analysis the maps are rotated to align the vertical axis (B) and then flipped to give a right-handed co-ordinate system (C). 



marginalized along <fi. This was done by discretizing 
both distance r and orientations a , f} and tf> and bin- 
ning the connections. Given the full distribution it was 
then possible to calculate the relevant marginal distribu- 
tions such as p(a \p, r) and p(fi\ a , r). 

Co-circularity 

If an injection site has orientation preference O i , then at 

a bouton site, located at angle <j> in the cortex relative to 
the injection point, there is a unique orientation prefer- 
ence Q£ c that is fully co-circular with the injection site 
(Figure 3B): 



9 b cc =2<t>-6 i (mod 180° 



(1) 



The difference in orientation between the actual 
orientation preference at the injection site 6 i , and the 

co-circular orientation Qy C , is a measure of the degree 

to which the bouton connection is co-circular. We 
define 



'diff 



I mod 18 



(2) 



We use the notation |-| mo d iso- to indicate that the 
absolute difference is calculated modulo 180°: 

I a | mod 180 .= min ( a ( mod 180° ) , ( -a ) ( mod 180° )) . (3) 



Geometrically, this means that when calculating the 
co-circularity of an edge we always consider the smallest 
angle between the edge and the fully co-circular orienta- 
tion 9 CC . It follows that edges are never separated in 

orientation by more than 90°(for instance, if 0 b = 135° 

and = 20 then ddia = 65) so the maximum possible 

value of (idiff is 90°. We calculated the mean d^ia 
between injection sites and bouton sites at distance r: 



D diff (r) 



' pairs 



'b Imod 180° 



N 



pairs 



E \9i + 9 b - 2<j> | mod 180 o 
pairs 



(4) 



(5) 



N 



pairs 



AiiffM has the intuitively appealing property that it 
is 45° when no co-circularity is present (that is, when 
the connections from the injection site have no sys- 
tematic dependence with orientation and angle) and 
decreases if the connections are connected co-circu- 
larly. To calculate D Aiii {r) we binned r in intervals of 
\Sk. The first bin, referred to hereafter as r = 0, con- 
tained pairs where 0 < r <0.75X, the r = \.Sk bin con- 
tained the pairs where 0.75X < r <2.25X and so on. 
This bin spacing balances the requirements of having a 
reasonable number of pairs in each bin to reduce 
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Figure 2 Orientation preference maps with connections. The orientation preference map for all cases with the injection site (white) and 
bouton traces (black) overlaid. The direction of the white injection sites indicates the preferred orientation at these points. The boutons indicate 
horizontal connections originating at the injection sites and connecting to the bouton sites. The orientation preference maps have been rotated, 
flipped and truncated (see methods) to create right-handed co-ordinate systems aligned with visual space and the overlaid boutons and 
injection sites have been similarly transformed after alignment. Scale bar is 1 mm. 




Figure 3 Definition of terms. (A) is a diagram of the various definitions of orientation and position between the injection site and the bouton 
site used in the methods. (B) shows a representation of eq. 1 showing co-circularity. Two edges are co-circular when they lie tangent to a 
common circle. Given an injection site with orientation 9,- there is a unique orientation preference 8^ c (for any <f>) that is exactly co-circular with 
the injection site. 

\ 
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noise, while allowing fine enough resolution to investi- 
gate trends with r. 

In addition, the distribution of co-circularity from 
each injection/bouton pair can provide useful informa- 
tion since various distributions of <i di ff, the co-circularity 
of individual bouton sites (see eq. 2), would result in the 
same mean -D diff (r). In addition to calculating the aver- 
age co-circularity £> diff (r), the values of d^is = | 9 { + O b - 

2^1 mod 180° for each pair at distance r were therefore 
binned to create a histogram H (r, d). The r bins were 
the same as used for calculating -D diff (r) and the d diii 
bins were spaced 15° apart. Since we were not interested 
in the absolute number of points at different distances, 
we calculated and plotted the conditional probability 
P(d\r). This allows a basic characterization of the distri- 
bution of co-circularity values. 

Co-circularity with offset 

A more general property than co-circularity is co-circu- 
larity with offset t. Co-circularity with offset r would 
arise if on a co-circular arrangement of edges, position 
was rotated independently of orientation. If lateral con- 
nections are co-circular due to co-circularity of natural 
scenes the maximal co-circularity strength should be 
found at r =0. To test this we also calculated co-circu- 
larity with offset using the definition: 

V . |0;+0 b -20-2r| modl8O o 
D m {r,r) = ^^ (6) 



N 



pairs 



for t = 0 to 80° in 10° increments. Since co-circularity 
is invariant under <j> — > (j) + 90° searching over this range 
covers all possible values for t . We also calculated 
77 best (r) for each value of r, which was defined as the 
value of r which minimized D Al ^r, r ). 

Testing for co-circularity 

In the original analysis of the tree shrew data it was 
demonstrated that lateral connections are denser along 
the axis of injection site orientation and connect prefer- 
entially to sites of similar orientation. This evidence 
demonstrates that lateral connections are likely to prefer 
co-linearity. In order to test whether the co-circularity 
we found in the horizontal connections was simply due 
to this co-linearity, we recalculated -D diff (r) (eq. 4), but 
excluded all points lying near the axis of the injection 
site. This was done by excluding all injection/bouton 
pairs with /J <30° or /J >150 . 



Controls and statistics 

It was important to have a reliable control to assess the 
significance of any co-circularity measured in the maps. 
We are interested in the significance of any deviation 
from 45° that occurs in -D diff (r). One comparison was 
simply examining the values of D dl [i(r, t ) at t * 0 
since these are expected by our hypothesis to be less co- 
circular then with the true origin (r =0). This means 

that we would expect r best to be near 0 when connec- 
tions are the most co-circular. 

We created 99 control cases {N = 99) for each animal 
by adding Gaussian noise (cr = 100 pixels) to the bouton 
positions and calculated .D di ff(r)y for each of these con- 
trols (/' denotes the index of the control). The added 
noise is high enough to overwhelm any co-circular pre- 
ferences in the connections since it is much larger than 
the wavelength of the maps, which means that it is lar- 
ger than the size of iso-orientation patches in the maps. 

We then used a permutation test to find if the value 
£> dif f(r) for each animal differed significantly from the 
control values. We tested -D diff (r) independently at each 
value of r. Permutation tests were used because we did 
not know what distribution the control values might be 
from and because we only had one true value to com- 
pare with the control distribution. We calculated the 
difference between the mean of the control values -D diff 
(r)j and the animal case D diS (r): 



T .-Ddiff (r) 
A Q( r ) = ^^ -D m (r). 



(7) 



N 



Our null hypothesis was that D dii ^r) was from the same 
distribution as the controls AtiffMy If the null hypothesis 
is true then the value of -D diff (r) can be exchanged with 
one of the controls D diS (r)j without affecting the expected 
value of AQ(r). There are N + 1 = 100 possible permuta- 
tions under these exchanges (since we also include the 
original permutation). We calculated AQ(r), , the value of 
AQ(r) for each of these permutations; using i as an index 
over all the permutations. We then calculated the signifi- 
cance of the true AQ(r) as the likelihood of this difference 
arising under the null hypothesis: 



,(r): 



^ =i +1 H(|AQ(r) j |-|AQ(r)|) 



(8) 



N + l 



where H (.) is the Heaviside function (H (x) = 0 for 
X < 0, H (x) = 1 for x > 0). We considered P <0.05 to 
be significant. 
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Natural scenes 

We used an existing image dataset [24] for all natural 
scene analysis. These images are uncompressed, 
deblurred and have linearized intensity values and thus 
should have no artifacts due to compression or other 
steps in the image acquisition. This image set was also 
used when co-circularity in natural scenes was originally 
demonstrated [6]. The images are 12-bit greyscale with 
a resolution of 1536 x 1024 pixels and an angular reso- 
lution of approximately 1 min of arc per pixel. The 
orientation and orientation strength of each pixel were 
calculated using steerable filters provided by matlabPyr- 
Tools [25]. After extracting orientation and orientation 
strength we calculated -D diff (r) for each image analo- 
gously to how it was calculated for lateral connections. 
The exception was that when averaging d di[{ we 
weighted by the orientation strength of the pixels in the 
pair and we excluded pixels with very low orientation 
strength (essentially un-orientated pixels). This was 
done differently because we did not have information 
on orientation selectivity for most the animal data. We 
then calculated a natural scene value by re-averaging the 
values of -DdiffM for each scene to create an average 
across all scenes. We did this double averaging to 
ensure that the calculation was not overly preferenced 
towards particular types of scenes which have more 
orientation density (such as tree scenes). 

Results 

Alignment and topography 

Figure 2 shows the aligned orientation preference maps 
with overlaid bouton and injection sites for all animals. 
In each case the vertical meridian was selected using pre- 
viously published results of tree shrew topography [21] 
and the map rotated to make it vertical. All pixels repre- 
senting receptive field positions in the ipsilateral visual 
field (that is, those to the right of the vertical meridian) 
were removed and these sections were not used in subse- 
quent analysis due to the very different magnification fac- 
tor in this region of the cortex. For analysis, the 
remaining region of VI was treated as isotropic with a 
constant magnification factor. We have previously 
demonstrated that this approximation does not signifi- 
cantly affect our measures of co-circularity [20]. 

There was a varying number of boutons sampled from 
each animal. Boutons positioned outside the VI orienta- 
tion map or on the excluded side of the vertical meri- 
dian were ignored. The number of included boutons for 
each animal is listed in Table 1. Although there was a 
large number of boutons for each case, each bouton 
cannot be considered an independent measurement as 
many boutons are clustered together near a single 
region of the map. 



Table 1 Number of included boutons for each animal 



Animal 


No. 


ts9509 


5411 


ts9514 


4974 


ts9529 


10708 


ts9531 


4821 



Co-circularity 

We quantified the co-circularity of the lateral connec- 
tions for each animal separately. Figure 4 shows the cal- 
culated D diti {r, t ) (see Methods) for each animal using 
the connections shown in Figure 2. D di{f is a measure of 
co-circularity, values below 45° indicate co-circularity. 
Each animal was treated separately because there were 
substantial differences between individuals. To establish 
whether variations in D difi were significant, each animal 
was compared with 100 control cases generated by add- 
ing Gaussian random noise (a = 100 pixels ~ 1.5 mm) 
to the bouton positions. The noise was much larger 
than the wavelength of the maps (approximately 35 pix- 
els), ensuring that boutons were moved significantly 
outside their original iso-orientation patch. Comparisons 
were made using nonparametric statistics (see 
Methods). 

Co-circularity with offset r is a generalization of 
co-circularity. Edges are co-circular when tangent to a 
common circle and co-circularity with offset t is when 
edges occur at an orientation r to the tangent of a 
common circle (so r = 0 is co-circularity). Co-circular- 
ity with non-zero r is not found in natural scenes and 
thus we would predict that if co-circularity in the lateral 
connections is significant, then it should be present 
more strongly than co-circularity with any non-zero off- 
set t . We measured co-circularity with offset at a range 
of non-zero offsets to test this hypothesis (Figure 4). 

These results indicate that the degree of co-circularity 
varies significantly between animals. However, previous 
work has demonstrated that lateral connections are elon- 
gated along their visual axis, indicating likely co-linearity 
[17]. Since co-circularity is a generalization of co-linearity, 
any co-linearity present contributes significantly to the 
degree of co-circularity. In order to clearly establish the pre- 
sence of co-circularity rather than just that of co-linearity in 
the lateral connection preferences, we re-calculated our 
results with all co-linear connections excluded. Figure 5 
shows the same measurements as the previous figure, but 
without co-linear connections. 

Case ts9509 is strongly anti-cocircular at longer ranges 
(Figure 4) and this tendency remains even if co-linear con- 
nections are removed (Figure 5). Case ts9514 is strongly 
co-circular, although when co-linearity is excluded this 
tendency becomes less strong (and it is anti-cocircular at 
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Figure 4 Quantification of co-circularity. Ddiff was calculated for all cases as a function of cortical distance, indicated in both millimeters and 
wavelengths X. The dotted black line marks D diff (r) = 45°, deviations below this line indicate co-circularity. Significance (P < 0.05) compared with 
the control cases is marked as '*' on figure. We also report Tbest M and the difference D di ff(r, 0) - D di ff{/, r b est M). (A) Case ts9509 had D di ff 
significantly different from the control case at r = 1.5, 3, 4.5, 6X (P <0.01). The values of r which minimized the D dlff at r = 1 .5, 3, 4.5, 6k were 10°, 
60°, 60°, 60° respectively, with differences D di ff(r, 0) - D dlf f(r, r best (r)) of 0.2°, 21.9°, 31.2°, 44.0° respectively. (B) Case ts9514 had D djff significantly 
different from the control case at r = 1 .5, 3, 4.5, 6X (P <0.01). The values of t which minimized the D djff at r = 1.5, 3, 4.5, 6X were 0°, 10°, 0°, 10° 
respectively, with differences D di ff(r, 0) - D di ff(r, tbest M) of 0°, 0.08°, 0°, 2° respectively. (C) Case ts9529 had D di ff significantly different from the 
control case at r = 1.5, 3, 4.5, 6X (P <0.01). The values of r which minimized the D dlff at r = 1.5, 3, 4.5, 6A. were 0°, 80°, 80°, 40° respectively, with 
differences D di ff(r, 0) - D di ff(r, Tbest M) of 0°, 0.7°, 2.9°, 9.3° respectively. (D) Case ts9531 had D d ,ff significantly different from the control case at r = 
1.5, 3, 6X (P <0.04). The values of r which minimized the D dlff at r = 1.5, 3, 4.5, 6X were 0°, 70°, 60°, 10° respectively, with differences D diff (r, 0) - 
Ddiffv", T"best 0) of 0°, 3.9°, 1 0.3°, 1 1 .5° respectively. 



large r). The two other cases show inconsistent results at 
different wavelengths, although case ts9531 shows strong 
co-circularity for the longest connections. 

Examining the results in the four cases at different 
values of t confirms that there is significant variation 
in the statistics of wiring between the four animals. Case 
9509 has values of r best (r) which are near 45°, while 



case 9514 shows r best {r) near 0° (the other two cases 
show inconclusive results). These results demonstrate 
that the variation between co-circularity and anti- 
cocircularity between the cases is significant and are an 
additional verification that the wiring on these connec- 
tions is not simply co-circular as our initial hypothesis 
predicted. As we discuss later, these findings indicate 
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Figure 5 Quantification of co-circularity excluding co-linearity. D^r) was calculated with all bouton sites which are co-linear with the 
injection site excluded (30° > /? < 150°). All cases where the true value was significantly different from controls were marked as '*'. (A) Case 
ts9509 had D di ff significantly different from the control case at r = 1.5, 3, 4.5, 6k (P <0.01). The values of r which minimized the Ddifr at r = 1.5, 3, 
4.5, 6k were 30°, 50°, 60°, 60° respectively, with differences D dlff (r, 0) - D diff (r, r best (/)) of 5°, 32.8°, 32.5°, 47.7° respectively. (B) Case ts9514 had D diff 
significantly different from the control case at r = 0, 1.5, 4.5, 6k (P <0.02). The values of t which minimized the D d iff at r = 1.5, 3, 4.5, 6k were 0°, 
20°, 20°, 30° respectively, with differences D diff {r, 0) - D dm {r, r best (r)) of 0°, 7.1°, 21.1°, 36.4° respectively. (C) Case ts9529 had D diff significantly 
different from the control case at r = 1.5, 3, 4.5^ {P <0.01). The values of t which minimized the D diff at r = 1.5, 3, 4.5, 6k were 40°, 50°, 70°, 60° 
respectively, with differences D diff (r, 0) - D diff (r, Tb es t W) of 12.4°, 9.5°, 3.1°, 19.1° respectively (D) Case ts9531 had D diff significantly different from 
the control case at r = 1 .5, 4.5, 6k (P <0.02). The values of r which minimized the D diff at r = 1.5, 3, 4.5, 6k were 10°, 60°, 60°, 10° respectively, 
with differences D dlff (r, 0) - D diff (r, r best (r)) of 2.0°, 3.9°, 16.0°, 1 1.5° respectively. 



that the results cannot be explained by minor variations 
in topographical alignment between the animals. 

In order to better understand the statistics of edge 
arrangements we also considered some other ways of 
examining edge statistics. Firstly, because multiple rfdiff 
distributions would result in the same mean value, we 
plotted the distributions of ddiff and compared the 



distributions with those found in natural scenes (Figure 
6). We found that in natural scenes the distribution of 
ddiff (averaged across about 4,000 natural scenes) has a 
small constant negative gradient. Although the slope of 
this distribution is small, this represents a sample across 
many scenes and within an small region of a scene stee- 
per gradients are often found. For the lateral 
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(A) ts9509 



(B) ts9514 





(C) ts9529 



(D) ts9531 




(E) natural scenes 




600 



800 



Figure 6 Distribution of rf diff values for the four cases compared with natural scene statistics (A-D) For each animal the probability 
distribution P (d^r) was calculated by binning into bins of 20° in d dlff and 1 5X in r. Because of the low number of connections the results have 
a significant amount of noise. (E) An analogous statistic was calculated for natural scenes. Although difficult to see at this scale (which matches 
that in panels A-D), there is small decline in probability with increasing d dif f so the mean value D diff (r) <45°. Because a large number of natural 
scenes could be counted and averaged there was very little noise in these results. 



connections the distributions appeared to be symmetric 
distributions centred (by definition) on the mean, D di{{ . 
This indicates that for each individual injection site the 
connectivity at a given wavelength tends to be centered 
on a particular value of co-circularity (that is, being either 
predominately co-circular or anti-cocircular or neutral). 

Additionally, we directly examined the probability 
distribution of edges in both natural scenes and 



lateral connections (data not shown). We found con- 
nections are isotropic for short-range connections (as 
also seen in the £> diff results). Long-range connections 
have strong biases for both orientation and direction. 
However, these statistics are probably due primarily 
to the previously studied elongation along the axis of 
the injection site orientation [17], rather than any 
specific arrangement of edges, as they are also present 
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in controls (where noise is added to the bouton 
positions). 

Discussion 

We examined lateral connections in VI and compared 
the statistics of their connections with the edge arrange- 
ment of natural scenes. These connections develop after 
eye-opening [2,9] and therefore may be driven by nat- 
ural scene input. As found previously, at short distances 
the lateral connectivity is isotropic (and probably mostly 
inhibitory [8]). However, long-range lateral connections 
are clearly anisotropic. We hypothesised that these con- 
nections might display similar connectivity statistics to 
edges in natural scenes. However, we found that the 
connectivity statistics showed a range of behaviour for 
the different individuals and both co-circularity and 
anti-cocircularity were present. This study is, as far as 
we know, the first to directly examine lateral connectiv- 
ity in VI for evidence of co-circular connections. 

Some limitations of our results are as follows. Firstly, 
due to archiving problems, the original data from the 
study of [17] was only available in fragmentary form and 
we were only able to reconstruct sufficient information 
for our analysis from four tree shrews. Secondly, the 
experimental data has a significant amount of noise due 
to error in alignment of optical imaging with histology. 
Because co-circularity is a second order property these 
errors are compounded when calculating co-circularity 
statistics. An additional possible source of noise is based 
on more theoretical considerations and inherent to any 
experimental study of these connections. Suppose the 
brain has some desired function for specifying lateral 
connectivity (for instance some spread of co-circularity), 
it is unlikely this function is realised without a signifi- 
cant level of noise in the connectivity. Because of this, 
even without any experimental uncertainty, it may be 
difficult to recover a good approximation to the original 
generating function. Assuming the noise in connectivity 
is independent between coloumns, this can be overcome 
by increasing the number of injection sites, however, in 
this study only four sites could be used. The restricted 
captive environment in which the tree shrews were 
raised was probably not as rich as a truly natural envir- 
onment, however, at the level of statistical characterisa- 
tion considered here these differences are not likely to 
be consequential. Because of these limitations, it is diffi- 
cult to make definitive conclusions from our results. 
However, this data is still able to provide some impor- 
tant insights into functional connectivity in VI. 

The variations between the four animals are not 
explainable by minor variations in topography. We 
examined the value of co-circularity for each of the ani- 
mals for various values of r , which is the value of co- 
circularity that would be present if the assigned 



topography was rotated by X . In order for, for instance, 
case ts9509 to be in agreement with case ts9514 they 
would have to contain a combined error of 45° in the 
assigned orientation of the topography. Other work has 
demonstrated that the topography of tree shrews has lit- 
tle inter-individual variability and consistent local struc- 
ture [26]. Additionally, our previous work has shown 
that the co-circularity statistics we used are robust to 
other deviations in topography such as anisotropy. 

Previous work [19] has shown that co-linearity alone 
does not fully explain the anisotropy in the original 
paper of Bosking et al. [17] and some broader form of 
connectivity is required. Here, we postulated that co-cir- 
cularity similar to that found in natural scenes would 
explain this. However, we found that there is a large 
variation in co-circularity statistics between different 
injections sites (and animals). In particular case ts9509 
shows strong anti-cocircularity, while case ts9514 shows 
strong co-circularity. 

Several other studies have examined horizontal connec- 
tions in VI in a variety of animals. Early work in tree 
shrew [27] found long-range (around 2.5 mm) lateral 
connectivity and some indications of orientation selectiv- 
ity in the connections, as did similar work in cats [28]. 
There is conflicting information on whether connections 
in VI are strongly iso-linear [17] (tree shrew) or mostly 
isotropic [8] (cat). However, there is significant evidence 
from psychophysics [11,12] that lateral connections play 
a role in some forms of co-linear facilitation. It has also 
been demonstrated that connections are refined during 
development [3] and there are correlations in co-linear 
neuron spiking [18]. In summary, there is strong evi- 
dence for co-linear lateral connections, but also indica- 
tions that other types of connections are present. 

Previous work has shown that an important functional 
role of lateral connections is inhibitory, such as iso- 
orientation suppression [29]. However, the majority of 
connections we have considered here, particularly at 
longer distances (r >500 ^m), are believed to be excita- 
tory [8,17]. Excitatory connections are known to be con- 
nected co-linearly and to be important for co-linear 
facilitation [11-13]. There is some evidence that VI cells 
in some species may also be involved in border-owner- 
ship computation [30], however, since we are not aware 
of any characterisation of the type of connectivity neces- 
sary for these computations, we did not test the proposi- 
tion here. 

Conclusions 

We suggest that there may be a role for both anti-cocir- 
cular and co-circular lateral connections. Edges in natural 
scenes are preferentially co-circular (with a small but 
consistent bias). For tasks such as co-linear facilitation 
[12], excitatory co-linear lateral connections are 
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necessary. This sort of facilitation is thought to be 
involved in contour integration [7]. However, from 
another point of view, it is precisely when input conflicts 
with what is expected, (in this case co-circularity), that it 
is providing novel and important information (indeed, it 
has been shown that high entropy areas of visual scenes 
are usually at the center of gaze [31]). From this point of 
view edges that are anti-cocircular are the most salient, 
and excitatory connections to facilitate their detection 
are a reasonable postulate. It is possible that the brain 
uses a combination of these strategies to both detect the 
unexpected while facilitating description of the more 
expected case (co-circularity). Our findings are consistent 
with such a bifurcated connection strategy. 

An interesting experimental followup to this work 
would be to examine the lateral connectivity of animals 
reared under unusual rearing conditions (as done for 
orientation maps in [20]). Since these connections 
develop after eye opening it is likely that these rearing 
conditions would have an effect on lateral connectivity 
(indeed previous work has indicated that strabismus 
during development has an effect on connectivity [3]). 
This could provide further insight into the develop- 
mental mechanisms and statistics underlying lateral 
connections and the degree to which these connections 
are influenced by natural visual scenes during 
development. 
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